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A METHOD FOR ESTIMATING THE COMPONENTS OF LIFT OF 
•WING -BODY COMBINATIONS AT SUPERSONIC SPEEDS 
By Warren A. Tucker 


SUMMARY 

The lift of wing-body combinations at supersonic speeds is obtained 
by means of an approximate linear ized-theory analysis which involves an 
assumption concerning the geometry of the body. Rectangular and tri- 
angular wing plan forms are treated; the wing and body may be at dif- 
ferent angles of attack. Also, the body may end at the wing trailing, 
edge or may extend farther back. The results are given in the form of 
simple generalized curves from which the lift of any specific configu- 
ration is readily obtained. 

A comparison with several groups of experimental data is presented 
to show that the method predicts, with acceptable accuracy, the lift on 
the wing and on the body. 

A question proposed by Ferrari concerning the optimum wing incidence 
for minimum total drag is investigated; it appears that, at moderate or 
high supersonic speeds, wing incidence is of little value in reducing 
the total drag. 

Several charts are presented to show how the several components of 
the total lift are distributed. 


INTRODUCTION 


In the past few years, the problem of the lift of wing-body con- 
figurations flying at supersonic speeds has been considered by several 
authors, who have approached the problem in different ways. Spreiter 
(ref. 1) has given the solution for low-aspect-ratio triangular wings 
mounted on slender bodies which end at the trailing edge of the wing. 
Ferrari (refs. 2 and 3) has developed an iteration method applicable to 
rectangular wings mounted on bodies which may extend behind the w'ing 
trailing edge. Ferrari ’s work is notable for being the' only one that 
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treats rigorously the lift on the part of the body behind the wing 
trailing edge (commonly called the afterbody). Although the method in 
principle requires iteration, apparently the first step in the calcula- 
tion procedure may give an answer accurate enough for engineering use 
(see refs. 3 and 4). Even with this simplification, however, and in 
spite of the fact that some shortening of the numerical calculations 
is possible (ref. 5), the work involved in the first step of the itera- 
tion process is still sufficient to make Ferrari's method one which is 
adapted better to the intensive investigation of the details of a 
particular configuration than to the examination of the over-all char- 
acteristics of a wide range of configurations. 

A more recent contribution has been made by Nielsen and Pitts 
(ref. 6). This work includes an exact solution, within the framework 
of the linear theory, for the pressure distribution due to unit angle 
of attack of the wing of a rectangular wing-body combination. The 
field of application of this method is similar to that of Ferrari's. 

An approximate method of wider scope has been proposed by Morikawa 
in reference 7- By considering limiting cases, and by making assumptions 
concerning the pressures on the body and on the wing, Morikawa finds the 
integrated lift on the wing and on the body for a wide range of param- 
eters. He presents results for wings of rectangular, delta, and clipped - 
delta plan forms but considers only those cases in which the body is at 
the same angle of attack as the wing and extends no farther back than 
the wing trailing edge (no afterbody). As Morikawa points out, such an 
analysis should have experimental verification before being too widely 
applied. 

The present contribution is also an approximate method but is 
based on different assumptions from those used by Morikawa. The 
principal assumption is that the body can be replaced by a flat plate 
in the plane of the wing; the lift on the wing and on the body of the 
resulting planar configuration are then readily calculated by standard 
methods. The method is applied to configurations having wings of either 
rectangular or triangular plan form, which may be at a different angle 
of attack from the body. The effect of the afterbody is also treated 
in an approximate manner, so that the body may end at the trailing edge 
of the wing (no afterbody) or may extend farther back. A correlation 
is made with available experimental data in order to assess the validity 
of the simplifying assumptions. 

0 

While the present work was being completed, another approximate 
method proposed by Nielsen and Kaatari (ref. 8) became available. In 
that paper, the lift on the wing is estimated by a modification of 
slender-body theory (ref. 1). In order to estimate the lift on the 
body, the body is imagined to be collapsed to a flat plate, as in the 
present paper, but only the lift on the body due to wing angle of attack 
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is considered; the lift on the body due to the wing in the presence of 
the body upwash is neglected. The effect of an afterbody is approxi- 
mated in the same general manner as in the present paper, with a 
slightly greater effective area being used. The method is applied to 
configurations having wings of rectangular, triangular, or trapezoidal 
plan form at the same angle of attack as the body. A comparison with 
a large number of experimental data is presented to show that the method 
predicts very closely the total lift of wing-body combinations. The 
method has also been used to predict the division of lift between wing 
and body, but the results are not yet available in published form. 


SYMBOLS 


a 

body radius 


s 

wing semispan 


c 

wing root chord 


m 

slope of wing leading edge ( 

see fig. 2) 


rectangular coordinates for 

field points 


rectangular coordinates for 

source points 

Se 

exposed wing area 

,1 

V 

q 

free-stream velocity 
free-stream dynamic pressure 

V, 

M 

free-stream Mach number 

% 

\ . _ \ 

(B = )fW 

°V 

wing angle of attack 

N 

V 

a B 

body angle of attack 


°u 

upwash angle 


% 

wing angle of incidence, ay 

- as 

w 

upwash velocity 



w B = a B V 



velocity potential 
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• <P 
C L W 


wing lift coefficient. 


Lift on wing 


c Lb 


body lift coefficient, 


Lift on body 
qS e 


C D 


drag coefficient. 


Drag 

q s e 


A aspect ratio of rectangular wing, — 

c . 

All angles are in radians, unless otherwise specified. 


ASSUMPTIONS AND METHOD 


At the outset, the original problem. was broken down into three 
separate parts, as proposed in reference 9 by Lagerstrom and Van Dyke 
(see fig. l). In figure 1, the lift of the complete configuration, 
with the wing at an angle of attack ay and the body at an angle of 
attack as, is shown as equal to the lift of the isolated body plus the 
lift of the two configurations on the right-hand side of the figure. 

In each of these two configurations the body is at zero angle of attack 
and extends to infinity ahead of the wing; in the one case the wing is 
at an angle of attack ay, and in the other case the angle of attack 
varies with the spanwise position along the wing, being equal to a^, 
the upwash-angle distribution around the isolated body at angle of 
attack aB- If the body is assumed to be approximately cylindrical . 
from the trailing edge of the wing to a point far enough ahead of the 
wing so that the influence of the body nose can be considered negligible, 
then in the notation of figure «2 the upwash angle has been shown by 
Beskin in reference 10 to be given by the following equation: 



The lift of the isolated body can be assumed known (ref. 11 or 12), so 
that the remaining problem is to find the lift on the wing and on the 
body of the two configurations at the right-hand side of figure 1. It 


is remarked that within the framework of the linearized theory the 
superposition scheme shown in figure 1 is not an approximation but is 
exact in the sense that the sum of the solutions to the three sub- 
problems is the complete solution and requires no iteration. The 
purpose of the scheme is to simplify the total problem by means of a 
breakdown into several basic problems. 

The remaining problem is now considered. A rigorous solution 
(except for Ferrari *8 iteration method and the work described in ref. 6) 
has so far not been published and when it is found it will almost 
certainly be of such a length that, like Ferrari *s solution, it will 
be best adapted to detailed analysis of a particular, configuration 
rather than to a broad examination of the field. Therefore, in order 
to obtain useful approximate values for the lift on the wing and on 
the body, two assumptions, are made. The first of these assumptions is 
t hat the cylindrical body can be replaced by a flat plate in the plane 
of the wing pjut ^ the upwash is still taken to be that around the cylin- 
drical body] . This assumption, of course, violates physical reality 
because the boundary condition on the true body surface is no longer 
met; the success of the approximation is best judged by a later com- 
parison with experimental results. The second assumption , which is 
necessary only when the body extends behind the trailing edge of the 
wing, is that the lift on the body is confined within the area shown 
in figure 3. This assumption is again only an approximation to the 
actuaT’s ituat ion , as can be seen by an examination of the data of refer- 
ence 13 , for example. 

The reason for making each of these approximations is the same; 
namely, to simplify the problem to the extent that useful approximate 
answers can be obtained without excessive calculation for fairly 
extensive ranges of the variables involved, in the present case, the 
lift on the body and on the wing was determined in the following manner. 
The wing-body combination was represented by a source distribution 
appropriate to each case (rectangular or triangular wing and wing angle 
of attack given by ay or a^) , and the velocity potential in each of 

the several pertinent areas, taking due account of the tip effect where 
present, was determined by the methods of Puckett and Eward (refs. 14 
and 15). The potential in each area was then evaluated at the down- 
stream limit of the area; this step resulted in the spanwise lift dis- 
tribution, which was then integrated (in most cases numerically) over 
the wing and over the body to give the net lift on each component. A 
typical case which illustrates the procedure is described in the appendix. 
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PRESENTATION OF RESULTS 


Calculations have been made, by the procedure described previously, 
of the lift on the wing and on the body for configurations having rec- 
tangular and triangular wings. The results were obtained in the form 
of generalized lift-curve slopes as functions of the generalized wing 
plan-form variable PA (for rectangular wings) or pm (for triangular 
wings), with the body-wing size ratio a/s as a parameter. For a 
particular value of a/s, the calculations were carried out only for 
values of PA or Pm greater than a certain minimum value determined 
by the position of the Mach line from the leading point of the wing- 
body juncture relative to the. body or to the wing tip (for the rec- 
tangular wing only). The limiting cases are shown in figure 4. 


The results of the calculations are presented in figures 5> 6, 
and 7 for the rectangular-wing case and in figures 8, 9 , and 10 for the 
triangular-wing case. The curves for the triangular-wing case have 
been extrapolated past the limit shown in figure 4 to Pm = 0, as indi- 
cated by the dashed parts of the curves. In most cases the method of 
extrapolation is obvious; in the case of figures 8(a) and 9(a) the 
fairing was guided by reasoning concerning the qualitative manner in 
which the lift is divided between the wing and the body for combinations 
of Pm and a/s outside the limit of figure 4(b). The fairing is best 
justified by the later comparison with experimental data, where for 
some points the faired parts of the curves were used. In each figure 
the following components are presented: 


PC LW 


a W 


55 * 

a B 

pc l b 


ay 


the lift on the wing due to wing angle of attack 


the lift on the wing due to body upwash 


the lift on the body due to wing angle of attack 


a B 


the lift on the body due to the wing in the presence of the 
body upwash 


The origin of these four components is evident from an inspection of 
the two configurations at the right-hand side of figure 1. The lift 
components are presented in this manner in order to preserve the 
generality of separate angles of attack for the wing and body. The 
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total lift for each part of the configuration is found by simply adding 
the components. Thus, if Clw(b) is defined as the total lift coeffi- 
cient acting on the wing in the presence of the body, then 


o 





and if 


CL Boo(W) 


is defined as the total lift coefficient acting on the 


semi- infinitely long body in the .presence of the wing, then 


' L Boo(W) 


if pc l b 



( 3 ) 


If the body is 
C Lb 0 > the lift 

the tota l lift 
then 


not semi-infinitely long, then to Cp B , ^ must be added 
of the isolated body. Thus, if Cp-p ^ is defined as 
acting on the finite body in th e presence of t he wing. 


0 



^ fr 




PCl b + 

a B 



( 4 ) 


Values for 


PC L . 


B r 


a B 


calculated for a body lift-curve slope of_2 based on 


body base area (refs. 11 and 12), are given in figure 11. This figure 
merely changes the reference area from the body base area to the exposed 
wing area, the reference area for all lift coefficients presented in 
t his naner . 


COMPARISON WITH EXPERIMENTAL DATA 


The experimental data available for comparison with the calculated 
results are by no means numerous but they are sufficient to permit some 
comparison. In view of the fact that the data from the various sources 
are not all presented in the same form and are not equally complete (in 
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some cases tests have been made at only one Mach number, in other cases 
not all the components of lift have been measured) several groups of 
data are discussed separately. Insofar as is feasible, each group of 
data is presented in its original form. 

Langley 9- inch supersonic-tunnel data .- Although this series of 
tests has not yet been completed, some of the results have been made 
available to the author in unpublished form. Briefly, the investiga- 
tion consists of tests at three Mach numbers of ten wings (three rec- 
tangular and seven triangular) and a single body, tested alone and in 
combination. In addition to measurements of total lift, the lift on 
the wing in the presence of the body is measured by an internal balance. 
For the triangular-wing case this is the only measurement that is yet 
available; measurements for the rectangular -wing case are complete. 

Each test was made at two values of Reynolds number, but the scale 
effect was small. The data presented herein are for the higher Reynolds 
number. 

The experimental and estimated results for the rectangular -wing 
case are compared in figures 12, 13, and 14. The lift coefficients 
are based on total wing area and the lift-curve slopes are per degree 
measure. The quant ities compared are identified in the (a) part of 
each figure. The agreement in most cases is acceptable. 

The comparison for the triangular -wing case is shown in figures 15 
and 16. The lift coefficients for this case are based on exposed wing 
area. The agreement is satisfactory. 

Sparrow missile data . - These data were obtained from tests of a 
13. 5-percent-scale model of the Sparrow 13-D (XAAM-N- 2 ) air-to-air 
missile at the Ordnance Aerophysics Laboratory at Daingerfield, Texas 
and have been presented in reference 16 . Both the wing and the tail 
of the Sparrow missile are triangular but have different apex angles. 

The test program, which covered a range of Mach numbers from 1 . 5 to 
2.5, included tests of the body-plus -wing and the body-plus -tail con- 
figurations and it is the data from these tests which are. compared with 
the present estimates. Forces were measured on the body alone, on each 
combination (body plus wing and body plus tail), and on the wing in 
the presence of the body.. For the body-plus -wing configuration, tests 
were made not only with the body at the same angle of attack as the 
wing but also with the body fixed at zero angle of attack while the 
wing angle of attack was varied. 

The appropriate comparisons are made in figures 17 and 18. All 
coefficients are based on th.e exposed area of the wings and the lift- 
curve slopes are per degree measure. In figure 17 the lift of the 
isolated body has been subtracted from the measured forces where 
appropriate. The agreement is, in general, satisfactory; the correlation 
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for the body-plus -tail configuration offers some hope that the method 
of estimation may work for configurations with no afterbody, and the 
results for the forces on the wing in the presence of the body show 
that the method fairly estimates the division of load. In addition, 
the acceptable agreement for the variable- and fixed-incidence cases 
indicates that the effect of the body upwash on the wing and on the 
body is correctly estimated. 

Other data . - The foregoing groups of data are the only ones avail- 
able that are sufficiently complete to permit a direct comparison of 
estimated and measured interference effects and components of lift. 
There are other data, however, for the total lift of bodies and wings 
in combination and the lift estimated by the present method can be 
compared with the measured values. Such a comparison, of course, is 
not so valuable as the preceding ones, but it does provide a check on 
the usefulness of the method for estimating over-all lift. The data 
presented in reference 8 and references l6 to 27 have been compared 
with the estimated values of total lift and the results are shown in 
figures 19 and 20. For consistency, all lift coefficients have been 
based on exposed wing area and the lift-curve slopes are per degree 
measure. The comparison is shown both including and excluding the lift 
of the isolated body, since in a strict sense the estimation of this 
quantity is not a part of the wing-body interference problem. For 
most of the cases the correlation is satisfactory, although there is a 
tendency for the estimated values to be low, particularly for the 
triangular -wing configurations. 


DRAG DUE TO LIFT 


An interesting problem was posed by Ferrari in reference 2; namely, 
what value of wing incidence will produce the minimum total drag for a 
given lift coefficient of the wing-body system? A related question, of 
course, is the following: By what percentage may the drag be reduced 

by using the optimum value of wing incidence? Ferrari worked out. the 
answer for a particular configuration in reference 2 but the numerical 
results are not applicable because of an error (later corrected in 
ref. 3) in one of the boundary conditions. 

The following results are obtained by following an analysis similar 
to that given by Ferrari in reference 2. The part of the total drag 
which is due to lift is given by the following equation: 


1 aC D 

Pc L 2 ' 


a l a I> + a 2 a B a W + a 3 a W 


( 5 ) 
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The optimum values of the wing and body angles of attack are given by 
the following equations: 


H/ 1 


2a^ac> - 


HCl 2a 4 2 


2^3,^ — ELgEL^ + EL^ 


3°l a 4 


and the wing angle of incidence is by definition 


X W ^W_ ^B_ 

P C L P C L P C L 


p ° l b„ ec Ll 


_o + 


0% 

+ 

as a w 


3Cl Bq + + 

a B a B 

°W a w 

B^L Bo | bc Lb • 

a B °iB a B 
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The preceding relations were vised with the data presented in fig- 
ures 6, 8 , and 9 to determine the optimum angles of incidence for two 
series of wing-body combinations. One series had rectangular wings and 
the other had triangular wings; both had long afterbodies. The results 
are presented in figure 21. The most obvious conclusion that might be 
drawn is that only at low values of 3A or 3m (which, for a given 
physical configuration, correspond to low -supersonic Mach numbers) and 
for bodies that are large relative to the wings (high a/s values) is 
the optimum incidence angle positive. However, the actual values of 
the optimum angles are very small. For example, at a Mach number of 
,2.0 and a lift coefficient of 0.1, a value of iw/3C;L equal to 0.1 
corresponds to an incidence angle of about 1°. Furthermore, calculations 
of the actual drag for the optimum configurations showed that the reduc- 
tion in drag was never greater than 1 percent of the drag due to lift 
for the corresponding zero-incidence configurations. The final con- 
clusion, therefore, is that at moderate or high values of supersonic 
Mach number wing incidence is of little value in reducing the total 
drag. 


DIVISION OF LIFT 


The manner in which the total lift is distributed is of prime 
interest in the wing-body problem. The data. presented in figures 5 
to 11 have been used to find the division of lift for wing-body com- 
binations with rectangular and triangular wings at zero incidence for 
two values of the parameters 3A and 3m, respectively, and the results 
are shown graphically in figures 22 to 25 . Certain general features 
are demonstrated by both the rectangular -wing and the triangular -wing 
configurations; namely, that relatively more of the total lift is 
carried by_the body for the higher values of 3A or 3m, that the 
interference lift on the body is relatively less for these higher 
values of 3A or 3m, and that the interference lift on the wing is 
primarily a function only of the relative body size a/s. 

The interference lift as a fraction of the total lift may also be 
examined without regard to the division between wing and body; this 
form of presentation is shown in figures 2 6 and 27. Again the 
rectangular -wing and triangular -wing configurations exhibit the same 
general characteristics: The interference effect is relatively greater 

for long-afterbody configurations and for low values of 3A or 3m, 
and the value of the relative body size a/s at which the interference 
effect is greatest is larger for low values of 3A or 3m. 
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CONCLUDING REMARKS 


A simple method is described for estimating the- lift at supersonic 
speeds of wing-body combinations having rectangular or triangular wings. 
The body may extend behind the wing trailing edge and the wing may be 
at a different angle of attack from the body. 

The results obtained are compared with experimental data; good 
agreement is shown. The results are then used to show that, at moderate 
or high supersonic speeds-, the use of wing incidence is of little help 
in reducing the total drag. The relative values of the different com- 
ponents of the total lift are shown in a series of figures. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX 

LIFT ON TRIANGULAR WING DUE TO BODY UPWASH FOR 0m < 1 


The first step in determining the lift on a triangular wing due to 
body upwash is to find the velocity potential cp at any point on the 
wing or (flat) body surface. The methods of references 14 and 15 are 
used and the setup for the necessary integration is shown in figure 28. 

The potential is given by the following equation: 


JTCp 

w B a 2 


^ m(x+0y) 

l+0m 


(a + q)‘ 


px+0(y-q) 


d! 


1_|3m (x+0y)+0q ^ 


(Al) 


I) 2 - 3 2 (y - I) 2 


l+0m 


which when integrated gives the following result: 


itcp 

w B a 


1 + 7 


cosh 


-1 


(f 


1 7) ♦ (f ♦ n~ f 7) 

|(1 + 0m)7 


+ 



(1 + + 3m | 7) f |(l + 0m) 


cosh - ’*' 


2(| + 0m J 7) + f(l + 0m) 
f(l + 3m) 


(A2) 


where 
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Now the spanwise lift distribution on the wing is given by the following 
expression: 


nV dL _ 

4qw B a 2 ^ ^ 


(A3) 


where is the velocity potential evaluated at the wing trailing 

edge. 


The total lift on one exposed wing panel is the following: 


L 

qa B 


4a 2 



Jo 


* 


$> a 7 


(A4) 


/ Wb \ 

^since — = a -qJ and the lift coefficient becomes 


or 


CLw _ 2L/ qa B 
a B S e 


8m 

it 




a/s 


0 d/ 


55 * 

a B 




0 d7 


(A5) 


(A 6) 


and when $ is evaluated from equation (A2 ) , then the final expression 
for the lift coefficient on the wing results: 
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P °LW 

a B 


1-* 

0 £ 


( a / s \ 2 

P-7* 

a/s 

^ cosh"^ 

2 ( 1 - 1) - ft 1 - Ml 

V-fj 


1 + 7 

|(1 + pm )7 

L 

0 

— 



^ - I - f 7 

(1 + 7)/l + 3m |(1 + 7) 


cosh 


J2 - |(1 - Pm) +■ 23m | 7 


a(i + 

S ' 


3m) 


! d7 (A7) 


This expression was evaluated numerically to give the curves shown in 
figure 8(b). 
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Figure 




(a) Lift on the wing due to wing angle of attack. 



(b) Lift on the wing due to body upwash. 

5.- Lift on wing for wing-body combinations with rectangular wings. 
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Figure 6.- 



(a) Lift on the body due to wing angle of attack. 



(b) Lift on the body due to body upwash. 


Lift on body for wing-body combinations with rectangular 
wings and no afterbodies. 
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(a) Lift on the body due to wing angle of attack. 



(b) Lift on the body due to body upwash. 

Figure 7 - - Lift on body for wing-body combinations with rectangular 

wings and long afterbodies. 
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Figure 11.- Lift of isolated body. Lift coefficient based on exposed 

wing area. 












NACA RM L52D22 


31 



Figure 12.- Comparison with Langley 9- inch supersonic-tunnel results 
for rectangular -wing case. Lift on wing due to "body. Coefficients 
based on total wing area; a in degrees. 
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(c) M = 2.kO. 

Figure 1 6 .- Comparison with Langley 9-inch supersonic-tunnel results 
for lift on triangular wing in presence of body. Sweepback angle 
of leading edge, 45°. Coefficients based on exposed wing 
“a in degrees. 


area; 




Measured 









Lme of perfect 
agreement 



Triangular wings. Coefficients based on exposed wing area; 
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(b) Triangular wings. 

Figure 21.- Angle of incidence for minimum drag due to lift. 
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